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Abstract 

We report the resuhs of a set of simulations of small arrays of submicron ferromagnetic particles. 

The actions of dipolar and exchange interactions were qualitatively investigated by analysing the 

O 

"^ ' ferromagnetic resonance spectra at 9.37 GHz resulting from the magnetization response of con- 

nected and unconnected particles in the array as a function of the applied dc magnetic field. We 
find that the magnetization precession movement (at resonance) observed in individual particles in 

Ch ' the array presents a distinctive behaviour (an "amplitude mismatch") in comparison to isolated, 

'^ ■ one-particle reference simulations, a result that we attribute to the action of interparticle dipolar 

O . couplings. Exchange interactions appear to have an important role in modifying the spectra of 

connected particles, even through a small contact surface. 
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OO ■ PACS numbers: 75.30.Ds, 76.50.+g, 73.21.-b 
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I. INTRODUCTION 

The main theoretical framework that defines micromagnetism has been estabhshed and 
refined several decades ago by Landau, Lifshitz and Gilbert, among othersii^i^. Micromag- 
netism describes magnetic phenomena at scales of ~ 10^^ — 10^^ m and offers a formal 
basis for the comprehension of emerging and collective phenomena in magnetic materials, 
as for instance magnetic domains in ferro magnets^, collective excitations of spins ("spin 
waves" )^i^>^, etc. It is based on a continuum, semi-classical dynamical model for the gyro- 
magnetic precession, including a phenomenological damping term. 

During several years^, micromagnetism has been somewhat overlooked due to difficul- 
ties in devising high-resolution (submicrometric) experiments to test it. In addition, the 
fundamental equation of micromagnetism, known as the Landau-Lifshitz-Gilbert equation 
(hereon, LLG)^, can only be solved analytically for special cases or supposing considerable 
simplifications for the description of the problen^. This situation changed rapidly in the 
past few years, with the development of improved experimental techniques and by the use 
of numerical simulations with exponentially increasing performance rates, allowing a deeper 
understanding of the physical phenomena encompassing the realm of magnetism at submi- 
croscopic scales^i^iii. These studies mainly aim at the application of microelectronics and 
data storageiS. In particular, the study of the role of the dynamics of confined spin waves 
in patterned arrays of magnetic elements in thin films^^'i^'i^'i^'i^ is of great interest in those 
applications, given that the excitations of spin waves effectively limit the time scale of the 
magnetization reversal process. 

The main purpose of this work is to study the ferromagnetic resonance (hereon, FMR) 
spectra emerging from small arrays of ferromagnetic particles, with different spatial config- 
urations and geometries, using as references single, isolated one-particle simulations. The 
idea is to identify qualitatively important contributions from dipolar and exchange couplings 
between spins as a function of the array configuration. In order to obtain an overall un- 
derstanding of these contributions, we have studied configurations varying from particles 
separated by a small distance, to particles slightly touching each other, or being arbitrarily 
connected with ferromagnetic material. The present simulations were shown to be com- 
putationally demanding, so at this time we have not surveyed the problem through an 
exhaustive number of array configurations, but instead focused on general trends from some 



representative cases. 

This paper is organized as follows. A theoretical background summary and the time 
domain micromagnetic simulation setups are given in Sec. I. In Sec. II, we describe the 
FMR spectra of the array configurations and the equilibrium magnetization fields. In Sec. 
IV, we discuss the results. 

II. MICROMAGNETIC SIMULATIONS 
A. Theoretical Background 

— * 

Let the magnetization vector M be defined as the sum of A^ individual magnetic mo- 
ments fij {j = 1, . . . , A^) in a small volume dV at a given position r of a ferromagnetic 
particle, namely: M{r) = — ^Jy ^ . Micromagnetism assumes that the direction of M varies 
continuously with position^. The dynamics of the magnetization field M{f,t) under the 
action of a external magnetic field Hext is that of a precession movement of M around an 
effective magnetic field Heff, defined as Heff = —f^o^ f."// ■ In this expression, Eefj rep- 
resents the energy associated with the effective magnetic field, and is given by the sum of 
four fields representing different interactions among the magnetic moments (spins) of the 
magnetic material, namely: E^fj = E^xch + Eanis + E^ag + Ezee- These four terms are, 
respectively: the exchange energy, the anisotropy energy, the magnetostatic or dipolar en- 
ergy and the Zeeman energy (namely, the energy associated with the external magnetic field 
Hext)- The equilibrium state of such a ferromagnetic system is that in which the total energy 
is minimized. 

The Landau-Lifshitz-Gilbert equation mathematically describes the dynamics outlined 
above, with the following additional physical parameters specified: the saturation magneti- 
zation, Ms (determined by the temperature, here fixed throughout), the gyromagnetic ratio, 
7, and a phenomenological damping constant, a. The LLG equation is then given by: 

^^^^^ = -^M{r, t) X H,jf{f, t) - ^Mir, t) x [Mir, t) x #,^^(f, t)] . (1) 

Note that the magnetization vector precesses around the H^ff field and looses energy to the 
environment in accordance with the damping constant, tending therefore to align with H^ff, 
given sufficient time span. 



The magnetization precession movement can proceed uniformly or not. The latter case 
leads to spin waves. Their occurence was confirmed experimentally and their behavior has 
been the object of several numerical investigations and comparative studies with the initial 
theoretical predictions^i^i^. We briefly summarize these possible magnetization precession 
movements in a small ferromagnet. The uniform precession movement of M{f,t) about the 
direction of an effective local field will occur at a given frequency ujq. The application of 
an oscillating magnetic field Hac at ujq, perpendicularly to the former field, will result in a 
coupling of M and Hac with energy absorption from the ac field by the system (leading to 
an "uniform resonance" or major peak in the FMR spectrum). But the ac field may also 
couple to nonuniform (spin wave) modes of precession of the M field. Exchange and dipolar 
interactions contribute to the energy of these modes. The former (exchange) interactions are 
expected to be more dominant at physically smaller magnetic elements, leading to additional 
resonances at frequencies Up = UQ+Dk^, according to the Kittel's model^, where D is a factor 
that depends on the exchange interaction between adjacent spins, and kp is the (quantized) 
wave vector corresponding to a given spin wave excitation in the ferromagnet. The resonant 
peaks associated with the exchange interactions lie at the left of the uniform peak— i^, 
that is, at smaller values of the dc component of the applied field. The later (dipolar) 
interactions, on the other hand, are relatively independent of the size of the ferromagnet, 
and may be important in a lattice of ferromagnetic elements^, leading to an interparticle 
dipolar coupling field. 

In order to probe the spin standing modes of coupled dots, two equivalent methods are 
available. The first one fixes the frequency of an small amplitude applied ac field for different 
values of a static magnetic field^. This is the method adopted in the present work (see next 
section). The second one fixes the static field and probe the frequencies of all the (quasi- 
uniform and non-uniform) modes. The former method is suitable to compare with FMR 
results, the latter with Brillouin scattering (BLS) results'^ . The behavior and character of 
the modes are the same. The latter method has been developed since about 2000 (pioneering 
work of Jorzick et al.^ ) and since then the dipolar-exchange nature and the symmetry of 
the spin modes of submicrometric dots has been fully understood. In particular, modes with 
nodal planes either parallel or perpendicular to the static applied field and edge modes have 
been identified, in addition to the quasi- uniform mode (see Gubbiotti et al.^^ and references 
therein). It is accepted that the modes with nodal planes parallel to the magnetization 



are high frequency modes (or equivalently on the "left" of the FMR peak) and modes with 
nodal planes perpendicular to the magnetization can exhibit frequencies lower and higher 
than the quasi-uniform mode (or, in other words, on the "right" and "left" of the FMR peak), 
depending on the number of nodes and the balance between the dipolar and exchange effects. 
The description of spin waves in a magnetic element can be given in analogy to a vibrating 
membranei^, but in the magnetic case two additional "restoring forces" take place, leading to 
a more complicated description of the normal modes than in the membrane case, where the 
description is made in terms of sinusoidal standing waves with a unique restoring force acting 
on the membrane. A thorough review of confined spin waves can be found in Demokritov 
et al.^. 

B. Simulations Setup 

We have performed micromagnetic simulations based on the numerical integration of the 
LLG equation (c.f. Eq. [1]) using the freely available integrator OOMMF (Object Oriented 
Micromagnetic Framework) ^^. In the present work, we have focused our investigations on 
a set of different small ferromagnetic (permalloy NisQFe2o) circular particle arrays with a 
small, finite thickness. Full (3D) simulations of the magnetization vector dynamics were 
performed, and interparticle magnetostatic or dipolar interactions were explicitly consid- 
ered in the computations. Simulations of isolated, single particles, were also performed for 
comparison purposes. 

The main global parameters of the OOMMF simulator are listed in Tab. [H and were 
fixed for all simulations here considered. Only parameters dependent on the specification 
of the particle geometries differ (these can be seen in Fig. [1]), as well as the number of 
external applied dc magnetic fields used to generate the FMR spectrum in each simulation. 
Table [IT] lists specific data of the simulations that may differ from each other, e.g.: the 
diameter of the particles (d), the diameter of the array (D), the interparticle spacing (a; 
i.e., the distance between the centers of two adjacent particles, or lattice spacing) and other 
simulation parameters, to be discussed in a moment. 

The arrays hold a small number of particles each, arranged in a regular, 2x2 (four 
particles, B- "family") or 3 x 3 (nine particles. A- "family"), grid (2D square lattice); see 
Fig. [TJ In the case of four particle arrays, we have studied three types of configurations: 



equally spaced particles (Bl configuration), equally spaced particles with arbitrary connec- 
tions among them (B2; in this case we define the interparticle spacing a as that of Bl, 
although the particles are connected among themselves)^, and particles "touching" each 
other (B3). As a reference, we have considered isolated particles (labelled ZO, AO and BO), 
which differ slightly in diameter, according to the particle diameters of the corresponding 
configurations, as indicated in Fig. [H Notice that we have generally prioritised the defin- 
tion of the array diameter with the choice of more "rounded" values for the discretization 
(defined by the "cell size" parameter; see Tab. [TTl) over individual particle diameters, but 
such a choice is immaterial; the effect of different relative particle diameters are considered 
in the analysis. The ZO particle simulation was included in order to compare with previous 
work by Jung et alM^. 

Larger arrays of particles were not considered at the present time, given the high com- 
putational demand of these simulations. Indeed, in order to obtain accurate results, the 



value of the cell size should not exceed the exchange length^, defined as lex = V'^^/ifJ'oMg), 
which in the present work results in /ex ~ 5.7 nm. Most of the simulations were executed on 
a 3 GHz Intel Pentium PC running Kurumin Linux and on a 2GHz Intel Core Duo running 
Mac OS X, and each simulation took between ~ 6 to 28 hours of CPU time, depending on 
the configuration. 

In order to obtain the ferromagnetic resonance of each of the configurations previously 
described, the following prescription was adopted. An external magnetic field in the plane 
of the particles was applied, formed by two components: a static (dc) magnetic field {Bdc = 
HoHfic) in the y direction, and a varying {ac) magnetic field {Bac = l^oHac) of small amplitude 
in the x direction, conforming with Jung et al.— : 

with the ac field frequency given by / = u;/(27r) = 9.37 GHz, A ~ /, and -Bac.o = 1 niT. 
Fig. [2] shows the time dependence and discretization of the applied Bac field. The values 
of the Bac field at intervals of 0.005 ns were used as inputs in the "field range" record of 
OOMMF, which in turn were stepped linearly by the simulator. The simulations were run 
up to 5 ns, resulting in 1000 outputs for each simulation (i.e., one simulation for each value 
of the Bdc field), as indicated in Fig. O We have performed a set of simulations for each 
array configuration by varying the range of the B^c field from 0.00 T to 0.39 T, at intervals 



of 0.01 T, resulting typically in 40 simulations for each configuration, as indicated in the last 
column of Tab. Ull The saturated regime is at B^c > 0.5 T and therefore outside our analysis. 
Some configurations present a few more than 40 simulations (BO, Bl); in these cases, the 
additional simulations were performed more fine-grainly around the resonance peak in order 
to evaluate whether significant divergence was found in the results (see next section). The Al 
configuration had a fewer number of simulations due to high computational demands; hence 
a more coarse-grained "sampling" of the underlying FMR spectrum was obtained in this 
case, as compared to the other configurations. Finally, we point out that the magnetization 
field of the particles was initially aligned to the same direction of the external B^c field. 

III. RESULTS 

In Fig. [3], we present the simulated time dependence of the spatially averaged magneti- 
zation vector M, in the x direction, normalized by the saturation magnetization (Mg). The 
various curves represent the resulting time dependence for each array /particle simulation 
evaluated at their respective FMR peak outputs, to be discussed below. The curves were 
offset for clarity. The AO simulation is not shown due to the fact of being very similar to 
the BO one, so it is omitted for clarity. The time dependence of the Bac field (arbitrarily 
normalized in the main figure and in the insets) is also included for comparison purposes. 
The insets in Fig. [3] correspond to zoom-in regions of the initial time evolution (left inset) 
and the steady state regime (right inset). Notice that all magnetization field responses are 
out of phase with the applied field after transient effects vanish, and the phase responses 
are practically identical in all cases, only differing in amplitude. Interestingly, the B3 con- 
figuration starts off at a different phase, but catchs up soon after the first cycle. From that 
figure it is already clear that the magnetization field of all array simulations have, on aver- 
age, a smaller response to the external field than the corresponding reference, one-particle 
simulations (represented by the BO simulation in the figure). 

In order to obtain the FMR spectra of the configurations, we proceeded as follows. The 
first 3 ns of all data have been excluded. For each simulation in a given configuration 
(i.e., for each applied dc magnetic field), the Fourier transform of the spatially averaged 
magnetization vector M, in the x direction, was obtained (as already mentioned. Fig. [3] 
refers to the results at the FMR peak only). The amplitude of the maximum Fourier peak 



at each B^c field was then obtained, resulting in the FMR spectra of Fig. HI In Fig. El 
we show the derivatives of the FMR spectra. Regarding the main body of Fig. HJ we 
have applied a spline fit to the reference single particle simulations' data (ZO, AO, BO), but 
maintained the individual data points of the array simulations. In the insets of Fig. HI 
however, all data has been spline fit to facilitate the comparison of the overall behavior of 
the curves. The derivative FMR spectra (Fig. [5]) were obtained by derivation (in intervals 
of 0.022 T) of the latter splined curves. 

We report the following three overall observations. First, the resonance uniform mode 
peaks of the reference single particle simulations show the expected trendM, namely, a shift 
in the peak position as a function of the diameter of the particle (the peak position shifts 
towards smaller values of the external field B^c as the particle diameter increases). The 
trend is very small, given that the particles differ only slightly in diameter. Second, (i) the 
resonant uniform mode peak position of the Bl array configuration is shifted towards smaller 
values of the external field B^c as compared to the reference (BO) one-particle simulation, 
even though the particles that compose the Bl array have the same diameter of the BO 
particle. Also, as already observed in Fig. [3l (ii) the amplitude of the uniform mode peak 
of the Bl configuration is smaller than that of BO. Finally, (in) the secondary peak at the 
right of the uniform mode peak in BO (also clear in ZO and AO, see arrows in Fig. [5]) is 
no longer apparent in Bl. These statements, namely, (i)-(m), are also applicable to the 
Al configuration, although the effects seem to be slightly more pronouced in this case, as 
one can notice by inspecting the derivative FMR spectra of Fig. [5] (see the dashed lines 
connecting the corresponding resonance peaks). Third, as the particles of the Bl array 
are arbitrarily connected (B2) or approach to the point of "touching" each other (B3), the 
resulting spectra significantly evolve away from the BO-Bl spectra: (a) the amplitude of the 
uniform mode peak of both B2 and B3 decrease and the overall spectrum gets more spread, 
with the presence of secondary contributions at the left of the uniform peak position (B2 
and B3), as well as to the right of it (in the case of B2). (b) The B2 uniform mode peak 
does not appear to shift in position relatively to the Bl peak, whereas the B3 does, towards 
higher values of the Bdc field. 

In order to better understand the observed characteristics of the FMR spectra, we anal- 
ysed in more detail the aspect of the magnetization fields at the uniform resonance peak. 
In Fig. [6] the "snapshots" of the magnetization vector field (at FMR) at four points of the 
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cycle {ut = 0, tt/2, tt, 3tt/2) are shown, being selected from a cycle around ~ 4 ns (when 
transient effects are over). Different pixel tonalities correspond to different values of the x 
component of the magnetization vector field, which in turn was subsampled to display an 
arrow for the average of 9 vectors per cell element. Notice that we have included a sinusoidal 
function at the top marking the four points corresponding to the selected "snapshots", so 
that immediately below each point of the reference curve the magnetization state at the 
cycle point can be directly observed and compared with that of the reference single particle 
simulation. 

The most obvious feature of Fig. [6] is that the individual particles composing the array 
configurations do not oscillate in synchrony with the corresponding single particle simula- 
tions. When averaging out the magnetization field over the array particles, the emerging 
response is lower than that of single particles, resulting in the observed lower amplitude 
of the uniform peak resonance of the spectra (Figs. H] and [5]). Indeed, one could natu- 
rally expect that the individual particles in the array configurations B2 and B3 would have 
their magnetization field evolving somewhat differently than that of the single particle sim- 
ulations, since these configurations physically join the particles of the array, removing the 
circular geometry characterization of the individual particles (forming effectively a "larger" 
single particle with a different geometry). But in the case of the Al and Bl configurations, 
the most suspicious agent causing the desynchronization would be an interparticle dipolar 
coupling field. 

In order to illustrate more clearly the behavior of the oscillating magnetization fields, 
we have synthetized the main information of the oscillating pattern into a "correspondence 
scheme". We have used the single particle simulations to associate a circular symbol of a 
diameter proportional to a given cycle point aspect (the pixel tonalities distribution) of the 
magnetization field. Fig. [7] illustrates the correspondence scheme adopted. We have used 
this correspondence scheme to recast the previous results (Fig. [6]) into the diagrammatical 
forms of Fig. [8] (excluding the B2 configuration, which cannot be represented in a simple 
manner in the scheme). Under this new representation, the underlying oscillatory patterns 
are more clearly displayed or synthetized. One can see that the suspected desynchronization 
is an "amplitude mismatch": the amplitude of the x-component of the magnetization as a 
function of position in the magnet, in each of the particles in a given array, is not similarly 
distributed among the particles. 



Another interesting point is the following. If one recast the above representation into 
a matrix representation (B configurations as 2 x 2 matrices and the Al configuration as a 
3x3 matrix), one can observe that the Bij elements behave in opposition to the B2J ones 
(j = 1,2). In the case of the Al matrix, one can notice a similar trend: the elements in the 
Alij and Al^j rows (j = 1, ..., 3) evolve in opposition to each other. The central row AI2J 
(j = 1, ..., 3) appears to evolve independently. 

IV. DISCUSSION 

The simulations of small arrays of ferromagnetic particles here performed indicate that 
their spatial configurations influence the resulting FMR spectra in distinctive ways. Here we 
discuss the possible effects of these spectra from the point of view two physical interactions: 
dipolar and exchange mode interactions. 

• Dipolar interactions between the particles: 

Al and Bl configurations: since in these cases the particles are not physically in 
direct contact, we may presume that the dipolar interaction is causing the moderate 
decrease in amplitude of the uniform mode peak in comparison to the one-particle 
reference simulations. Indeed, one can suppose that the finiteness and symmetry of 
the arrays over which the particles are distributed favor an interparticle dipolar inter- 
action for which the combined effect is that of a mismatching of the local amplitudes 
attained by the magnetization field precession movement (observed at the uniform 
mode outputs). Another observation is the shift of the uniform mode peak towards a 
smaller value of the applied dc magnetic field, in comparison to the one-particle ref- 
erence simulations. It mimmics the expected spectrum of an effective single particle 
of larger diameter. One may notice that in the case of the Al configuration the whole 
spectrum seems to shift accordingly (see dashed lines of Fig. [5]); this does not seem to 
be the case of the whole Bl spectrum. The effect, clearly seen in the Al configuration, 
has already been experimentally measured with BLS and calculated by Gubbiotti et 
al.^ with the second alternative method mentioned in Sec. II. A for a similar 3x3 
array. The shift towards lower values of the static field of the modes corresponds to 
the increase of the mode frequencies in Fig. 3 of Gubbiotti et al. A direct compar- 
ison between the modes found in the simulation of Gubbiotti et al. and the ones in 
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the present work is difficult because in the former case the dynamic magnetization is 
shown, while in this work (c.f. Fig. [6]) the total magnetization (static plus dynamic) 
is shown. We intend to proceed these investigations in a future work. Regarding the 
nonuniform modes of precession, one can observe the following facts. First, the lowest 
energy mode (at the right of the uniform mode), seen in the reference-one particle sim- 
ulations (see arrows in Fig. [5] over the derivative spectra of ZO, AO and BO), diminish 
for the Al and Bl configurations (one cannot affirm that they completely vanish, given 
the smothing over a limited resolution in B^c used to derive the spectra). Second, a 
relatively larger gain in response of the nonuniform modes at the left of the uniform 
mode peak is found (c. f. main body of FigE]). This will be discussed in more detail 
below, since these peaks are thought to arise from exchange interactions, although will 
possible contributions from dipole interactions. 

B2 and B3 configurations: In these cases, the particles are directly in contact. In 
the case of B3 configuration, the "amplitude mismatch" effect seen in the previous 
cases is again clearly observed (see Figs. [6] and [8]). Yet, from these figures one can 
qualitatively observe that, looking at individual particles in the array, the magneti- 
zation field precess in smaller amplitude as compared to the BO, Bl configurations, 
although in similar pattern as the latter, resulting in a smaller relative amplitudes of 
the average magnetization field, which one can conffim from an inspection of Fig. [31 
Although we have not traced an ocillatory pattern scheme for the B2 configuration, it 
is remarkable that it somewhat seems to follow the pattern observed in the Bl config- 
uration (as one can observe looking only at the central parts of the particles of the B2 
array, excluding the dynamics occuring at their connections, c.f. Fig. [6]), from which 
the latter differs from the former only by the application of arbitrary connections. 
The amplitude of the average magnetization field of the B2 configuration is greater 
than that of the B3 configuration, possibly because the contribution coming from the 
central region of the particles of the B2 configuration can follow somewhat that of the 
Bl ones, yet being significatively smaller than the latter ones because of the averaging 
over the field behavior of the connected parts. Evidently, the resulting spectrum will 
be more complex in this case, with, for instance, the presence of a new nonuniform 
peak at the far right of the uniform one (see the arrow over the B2 curve of Fig. [5]). 
Another distinct aspect of the B3 configuration is the shift of the uniform mode peak 
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towards a larger value of the applied dc magnetic field, in comparison to the one par- 
ticle reference simulation, contrary to the trends of the Al and Bl, B2 configurations, 
in which the peaks move to the opposite direction. The origin of the B3 shift is also 
unclear at this point and deserves further investigation. 

• Exchange mode interactions: 

Al and Bl configurations: In these cases, the first nonuniform peak at the left of 
the uniform one does not appear to suffer significant changes in amplitude as com- 
pared to the reference one-particle simulations. However, as already mentioned, the 
spectrum of the Al configuration is clearly shifted as a whole to the left. The other 
higher harmonics to the left of the uniform mode from both Al and Bl seem to have 
slightly higher amplitudes in comparison to the one-particle reference simulation (spec- 
tra between B^c = 0, ...,0.08 T). Since the particles in these configurations are not in 
physical contact, the energy of these nonuniform peaks, understood to be mainly the 
result of exchange mode interactions, must have been suplemented by the injection 
coming from other means. In this case, dipolar interaction is a reasonable candidate, 
since the particles do not touch each other. This assertion would agree with the hy- 
pothesis raised by Jung et al.^^ on the possible coupling of the dipolar interactions 
with exchange spin wave modes. 

B2 and B3 configurations: In these cases, as expected, the spectra evolve in a 
complex way in comparison to the reference simulations. In particular, one can observe 
a relatively large gain in response of the nonuniform modes (Fig [5]). As observed 
previously, a fraction of this gain should be a result of the coupling of the dipolar and 
exchange interactions. But since the particles of these arrays are in physical contact, 
the exchange interactions should have a more important role. One remarkable aspect 
is the form of the B3 configuration spectrum which is quite different from that of 
the Bl configuration, and in some aspects approach more closely to that of the B2 
configuration (c.f. Fig. [5]). This fact implies that propagative effects arising from 
exchange interactions, even when are able to act through a small contact surface, can 
result in significant modifications in the FMR spectrum. 

Finally, an important remark should be mentioned. According to the work of Gubbiotti 
et al.^ , the effect of the interparticle dipolar coupling is to split the modes, spreading them 

12 



into bands, in the limit of large arrays. In the case of a 3 x 3 array, each mode should give 
rise to 9 modes, including degeneracy. Gubbiotti et al. have shown that this broadening is 
appreciable, even large, for the quasi-uniform mode. However, this effect is not visible in 
the present study (c.f. cases AO-Al and BO-Bl; Figs. |l]and[S]). The reasons for the lack of 
splitting are currently under investigation. 

In summary, in the present paper we have presented a set of 3D simulations of small 
arrays of ferromagnetic particles supposed here to represent small isolated sections of a 
patterned thin film. We have analysed the resulting FMR spectra and the magnetization 
field behavior at the resonance modes. We show that the spatial configurations and ge- 
ometries of the particles in the arrays influence the resulting FMR spectra in distinctive 
and perhaps unanticipated ways. We have attempted to isolate the action of dipolar and 
exchange interactions by studying arrays with particles both connected and not connected 
among themselves. These interactions appear to have an interesting role on the dynamics of 
the magnetization precession among the particles in the array, as described in detail in this 
paper (synthetized in Fig. [8]). Other simulations are intended to be performed in a future 
work. In particular, it would be interesting to perform simulations with the applied field at 
different nominal angles (specially, at 45°), using the same arrays here analysed. According 
to measurements of Jung et alM, the typical low energy peak of one-particle simulations 
appear more pronounced in these cases. Since this low energy peak is suspected to arise 
mainly from dipolar interactions, it would serve as an interesting tracer of the role of these 
interactions in the array. 
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TABLE I: Main global parameters adopted for the OOMMF simulator, fixed for all simulations 
in the present work, except for the "particle width/height" and "cell size" parameters, which are 
listed separately in Tab. II, as indicated. 

Simulation Parameter/Option Parameter Value/Option 



Saturation magnetization [A/m] 
Exchange stiffness [J/m] 
Anisotropy constant [J/m^] 
Damping constant 
Gyromagnetic ratio [m/(A.s)] 
Particle thickness [nm] 
Particle width/height 
Cell size 
Demagnetization algorithm type 



8.0 X 10^ 
1.3 X 10-11 
0.0 

2.21 X 10^ 

85.0 

see D in Table II 

see Table II 

magnetization constant in each cell 



"The value of the damping constant here adopted is far larger than the real one for Permalloy (< 0.01). 
A small value of damping constant would allow a better resolution of the absorption lines within the FMR 
spectra but would lead to prohibitive computation times. 
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TABLE II: Particular parameter values of the simulations: the diameter of the particles (d), the 
diameter of the array (D; for individual particles, D = d.), the interparticle spacing (a, when 
applicable), and cell size. Other data are: the number of ferromagnetic particles (column 2) and 
the number of external applied dc magnetic fields (or number of simulations) used to generate the 
FMR spectrum (column 7). 

Label # FM particles d {fim) D (fim) a (pm) Cell size (nm) B^c (number of simulations) 

40 
40 
46 
26 
43 
40 
40 



zo 


1 


0.500 


0.500 


- 


5.00 


AG 


1 


0.559 


0.559 


- 


5.59 


BO 


1 


0.591 


0.591 


- 


5.91 


Al 


9 


0.559 


1.900 


0.671 


5.00 


Bl 


4 


0.591 


1.300 


0.709 


5.00 


B2 


4 


0.591 


1.300 
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FIG. 1: Particle configuration geometries, arranged in "families" of common characteristics. 
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FIG. 2: The time dependence and discretization of the applied Bac field. 
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FIG. 3: Time dependence of the amplitude (at FMR) of the spatially averaged magnetization 
vector Mx/Ms of the array configurations. In the main body of the figure, the curves were offset 
for clarity, and are shown at the same presentation order as the legend (top curve is that of the Al 
configuration). The time dependence of the external Bac held (arbitrarily normalized; thin line) is 
included. 
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FIG. 4: The ferromagnetic spectra of all configurations. 



19 



to 



ro 



to 




0.1 0.2 0.3 

B_dc [T] 



0.4 J 



FIG. 5: The derivative ferromagnetic spectra of all configurations. 
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FIG. 6: "Snapshots" of the magnetization vector field (at FMR) at four points of the cycle (wt = 
0, vr/2, vr, 37r/2; see sinusoidal curve at the top of the snapshots) after ~ 4 ns (when transient effects 
are over). Different pixel tonalities correspond to the value of the x component of the magnetization 
vector field at the pixel element. 21 
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FIG. 7: The oscillatory pattern correspondence scheme (see the main text form further details). 
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FIG. 8: Same as Fig. 6, recasted according to the scheme of Fig. 7, synthetizing the combined 
effect of the local amplitude of the x-component of the magnetization field at each cycle point of 
the ac driving field, for configurations Bl, B3 and Al. 
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^^ Our motivation for considering the B2 configuration was the study of the relative contribution of 
dipolar and exchange interactions for complex geometries. In addition, as one can observe in the 
paper of Jung et al.— (their Fig. 1), the patterning procedure appears not to be perfect and often 
produces artifacts between the dots. We have attempted to include arbitrary "bridges" between 
elements, with no particular criteria for their geometrical shape. Finally, the understanding 
of the micromagnetics of such complex geometries is being currently conducted for specific 
applications of interest at the AMR/IAE/CTA Division. 
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